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Abstract. We extend the formalism of weak interaction processes, obtaining new expressions 
Q_) , for the transition rates, which greatly facilitate numerical calculations, both for neutrino-nucleus 

^ ■ reactions and muon capture. We have done a thorough study of exclusive (ground state) 

C ' properties of ^^B and ^^N within the projected quasiparticle random phase approximation 

(PQRPA). Good agreement with experimental data is achieved in this way. The inclusive 
ncutrino/antineutrino (I'/t') reactions '^^C(!^, e~)^'^N and ^■^C(i',e^)^^B are calculated within 
^ I both the PQRPA, and the relativistic QRPA (RQRPA). It is found that the magnitudes of the 

QQ . resulting cross-sections: i) are close to the sum-rule limit at low energy, but significantly smaller 

than this limit at high energies both for v and i>, ii) they steadily increase when the size of the 
C^ , configuration space is augmented, and particulary for v/iy energies > 200 MeV, and iii) converge 

CO ■ for sufficiently large configuration space and final state spin. 
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1. Introduction 

The neutrino-nucleus scattering on ^^C is important because this nucleus is a component of 
many liquid scintillator detectors. As such it has been employed in experimental facilities 
C^ . LSND, KARMEN, and LAMPF to search for neutrino oscillations, and for measuring neutrino- 

nucleus cross sections. The ^^C target will be used as well in several planned experiments, 
such as the spallation neutron source (SNS) at Oak Ridge National Laboratory, and the Large 
Volume Detector (LVD) at Gran Sasso National Laboratories. On the other hand, this nucleus 
is important for astrophysics studies, as it forms one of the onion-like shells of large stars 
before they collapse. Concomitantly, the LVD group have stressed recently the importance of 
measuring supernova neutrino oscillations. A comprehensive theoretical review on neutrino- 
nucleus interaction was done by Kolbe et al. [1]. 

2. Formalism for the Weak Interacting Processes 

The most widely used formalism for neutrino-nucleus scattering was developed by the Walecka 
group [2]. They classify the nuclear transition moments as Coulomb, longitudinal, transverse 
electric, and transverse magnetic. This terminology is not used in nuclear /3-decay and ^- 
capture, where one only speaks on vector and axial matrix elements with different degrees of 
forbiddenness: allowed (GT and Fermi), first forbidden, second forbidden, etc. Motivated by 
this fact, we extend the formalism of weak interaction processes developed in [3, 4], obtaining 
new expressions for the transition rates, which greatly facilitate numerical calculations, both for 
neutrino-nucleus reactions and muon capture. 



The weak Hamiltonian is expressed in the form [2] 

H,4r) = -^J„^,e-^^■^ (1) 

where G = (3.04545 it 0.00006) x 10^^^ is the Fermi couphng constant (in natural units), the 
leptonic current /„ = {1, i/0} is given by [4, Eq. (2.3)] and the hadronic current operator 
Jq, = {J,jj0} in its nonrelativistic form [4]. The quantity k = Pi — Pj = {k, ^^0} is the 
momentum transfer, M is the nucleon mass, and Pi and Pf are momenta of the initial and 
final nucleon (nucleus). The effective vector, axial- vector, weak- magnetism and pseudoscalar 
dimensionless coupling constants are defined in [4]. In performing the multipole expansion of 
the nuclear operators Oa = (O,iO0) = J^e"*'''^ it is convenient: 1) to take k along the z axis, 
i.e., e~ —7- e"***^, and 2) to make use of Racah's algebra. One obtains 

O0J = gvM] +igAM^ + iig^ + gp,)M^j (2) 

Om^ = i{Smogp2 - gA + mg^)M:^j + gvMl^j - Smog^M], (3) 

with the elementary operators 

M] = ij(p)lj(r). Mi, = J2 i'-^-'Fumjiip) [Yar) €D ct]j , 

L>0 

M^ = M-'j,{p)Uf){a . V), M^^^j = M-i ^ i'-^~' FunJiipWar) V]j, 

L>0 

(4) 

where -Fum are Clebsch-Gordan coefficients [4]. The CVC relates the vector-current pieces 
of the operator Oa as k • O^ = kOq = kf^O^ with kg) = k^^ — 5(A£^coui — ^M), where 

AEcou\ — ^§jf- — 1.45Zyl^^/^ MeV, is the Coulomb energy difference between the initial and 
final nuclei, AM = M„ — Mp = 1.29 MeV is the neutron-proton mass difference, and S = ±1 for 
neutrino and antineutrino scattering, respectively. 

The transition amplitude 7jjj(k) = ^^ ^^ |(J^|i7^y(/i;)|0^)| for the neutrino-nucleus reaction 
at a fixed value of n, from the initial state jO"*") in the {Z, N) nucleus to the n-th final state | J^) 
in the nucleus {Z ±1, N ^1), reads 

rj.(K) = AnG^i Yl KJ^||O«j(^)||O+)|2£,-2K((JJ^||O0j(^)||O+)(J-||Ooj(^)||O+)*)£0o], 

a=0,O,±l 

(5) 

where the momentum transfer is A; = p^ — q^, with p£ = {p£,iE£} and q^, = {c[u,iEu}. The 
lepton traces C^j, Cq, C±i and £00 are defined in [4]. 

3. Neutrino (antineutrino)-nucleus cross section 

The exclusive cross-section (ECS) for the state |J^), as a function of the incident neutrino energy 

Eu, is 



I I ZT* /"l 

ae{Jl,E,) = \^-lF{Z + S,Ee)J_ d{cos e)Tj^jK) 



(6) 



where E^ = E^ — wjtt and jp^j are the energy and modulus of linear momentum of the lepton 
£ = e,p, and ojjn = —k^n = Ey — E^ is the excitation energy of the state |JJ^) relative to the 



state |0^). Moreover, F{Z + S,Ei) is the Fermi function for neutrino {S = 1), and antineutrino 
{S = —1) processes, respectively. Here, we will also deal with the inclusive cross-sections (ICS), 






as well as with folded cross-sections, both exclusive, and inclusive 

a^J^) = j dE^aei^l, E^)ni{E^) , ^i = J dE^ae{E^)ne{E^), 



(7) 



(8) 



where ni{Ey) is the neutrino (antineutrino) normalized flux. In the evaluation of both neutrino, 
and antineutrino ICS the summation in (7) goes over all n states with spin and parity S'^ < 7^ 
in the PQRPA (Projected Quasiparticle Random Phase Approximation), and over ^^ < 14 in 
the RQRPA (Relativistic QRPA). 

4. Numerical results 
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Figure 1. (Color online) ECS ae-{lt,Ey) for the reaction ^'^C{i^e,e^)^'^^ (in units of 10"^^ 
cm^), as a function of the incident neutrino energy E^. On the left side we have t = for all 
5jv, whereas for the right side i = for 5*2, and 5*3, t = 0.2 for ^4, and t = 0.3 for 5*6. The 
experimental data in the DAR region are from Ref. [6]. 



The calculations were performed with the PQRPA using the model spaces 52,^3, ^4, and 
Sq including 2, 3, 4, and 6 hcj harmonic oscillators shells, respectively, and by employing a 
(5-inter action. For 52, ^3, and 54 spaces the s.p. energies and pairing strengths were varied in a 
X^ search to account for the experimental spectra of odd-mass nuclei ^^C, ^^B, ^^C, and ^^N [4]. 
As this method can not be used for the 5*6 space, which comprises 21 s.p. levels, the energies 
in this case were derived as in Ref. [10], while the pairing strengths were adjusted to reproduce 
the experimental gaps in ^■^C [11]. We also employ the RQRPA theoretical framework [5] with 
520, and ^20 spaces. In this case, the ground state is calculated in the Relativistic Hartree- 
Bogoliubov model using effective Lagrangians with density dependent meson-nucleon couplings 
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Figure 2. (Color online) Inclusive cross sections a^^-^Ey), and ag^+{Ey) (in units of 10~^^ 
cm^) for the reactions ^^C(fe, e~)^^N, and ^^C(t'e) e"'")^^B, respectively, plotted as a function of 
incident neutrino and antineutrino energies. The PQRPA results within the s.p. spaces 5*2, 5*3, 
and 5*6, have the same parametrization as those on the right panel of Fig. 1. The sum rule limit 
SRgt and several previous RPA-like calculations, namely: RPA [8], CRPA [9], and RQRPA 
within ^20 for two-quasipartile cutoff E2qp=W0 MeV [10], as well as the shell model result [8], 
are also exhibited. 



and DD-ME2 parameterization, and pairing correlations are described by the finite range Gogny 
force. 

Fig. 1 shows ECS ^'^CgsiveT^^Y'^^gs plotted as a function of the incident neutrino energy 
El, for different final energy region E^ < 60 MeV (DAR region), E^ < 100 MeV ( supernovae 
neutrino signal search) and Ei, < 250 MeV ( neutrinos oscillation search in LSND). On the left 
side the ^-interaction particle-particle strength is t = for all Sjy ^, whereas on the right side t 
is gauged to reproduce the ground state energy of ^^N and the B{GT) values of ^^N and ^^B, 
getting t = for S2 and ^3, t = 0.2 for S'4, and t = 0.3 for Sq. The ECS experimental data in 
the DAR region are from Ref. [6]. 

On the left and right panels of Fig. 2 are displayed the ^^C(z/e, e~)^^N and ^^C(t'e, e"'")^^B ICS 
(Te=F(-E'iy) , respectively. The PQRPA results within the s.p. spaces 52, ^3 and S^, have the same 
values of t as in the right panel of Fig. 1. They are compared with: i) the sum rule limit SRgt 
obtained with an average excitation energy oTpr of 17.34 MeV, ii) other RPA-like calculations, 
and iii) the shell model (SM). 

In Fig. 3 are plotted the ag^{Ey), evaluated in RQRPA with different configuration spaces 
according the values of quasiparticle energy cut-off Egp. One sees that the ICS increase as the 
configuration space increase, saturating at higher and higher energies. 

Finally, Fig. 4 presents the RQRPA results for agT (-£■!/) with different maximal spins J. One 
sees that while for antineutrinos contribute spins only up to J = 11, it is needed to go up to 
J = 14 in the case of neutrinos. It is also clear that to achieve the saturation one should go to 
energies larger than 600 MeV. 



This value for 5*3 was adopted in Ref. [4]. 
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Figure 3. (Color online) Inclusive cross sections "^^C(fe, e~)"^^N (left panel) and ^■^C(t'e, e^)"^^B 
(right panel) (in units of 10~^^ cni'^) evaluated within the RQRPA for different configuration 
spaces. 
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Figure 4. (Color online) Inclusive cross sections ^^0(^6, e~)"^^N (left panel) and "'^^C(Pe, e"^)^^B 
(right panel) (in units of 10~^^ cm^) evaluated within the RQRPA for E2qp = 500 MeV and 
spaces S'20, and 530, for different maximal spins, which go up to J = 14 for neutrinos, and up 
to J = 11 for antineutrinos. 



5. Summary 

The shell model and the PQRPA are proper theoretical frameworks to describe the ground state 
properties of ^^B and ^^N and the ECS a^-{lt,E^) for the reaction ^'^C{iye,e-)^'^N. The ICS 
calculated, at difference with the exclusive ones, steadily increase, and particulary for neutrino 
energies larger than 200 MeV, when the size of the configuration space is augmented, in spite of 
including the particle-particle interaction. They approach to those of the first-forbidden sum- 
rule limit at low energy, but are significantly smaller at high energies both for neutrino and 
antineutrino. The study of the partial ICS's has been related with the proposal done in Ref. [7] 
of performing nuclear structure studies of forbidden processes by using low energy neutrino and 
antineutrino beams. We conclude that to study forbidden reactions in ^'^C{v',e~^y^B process, 
one would need P-fluxes with £";> up to > 150 MeV, while those from Ref. [7] go up to 80 MeV 
only due to the feasibility of the proposed experiments. 
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